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ABSTRACT 

We examine the effects of cool bottom processing (CBP) on the isotopic ratios 
18 o/16o^ 17o/16o^ C/0, and N/0 in the convective envelope dur¬ 

ing the thermally-pulsing asymptotic giant branch (TP-AGB) phase of evolution in a 
1.5Mq initial-mass star of solar initial composition. We use a parametric model which 
treats extra mixing by introducing mass flow between the convective envelope and the 
underlying radiative zone. The parameters of this model are the mass circulation rate 
(M) and the maximum temperature (Tp) experienced by the circulating material. The 
effects of nuclear reactions in the flowing matter were calculated using a set of static 
structures of the radiative zone selected from particular times in a complete stellar 
evolution calculation. The compositions of the flowing material were obtained and the 
resulting changes in the envelope determined. No major shifts in the star’s energy bud¬ 
get occur from the imposed CBP if log Tp < 7.73. Using structures from several times 
on the TP-AGB, it was found that the results for all species except ^®A1 were essen¬ 
tially independent of the time chosen if logTp > 7.6. Abundant ^®A1 was produced 
by CBP for logTp > 7.65. While ^®A1/^^A1 depends on Tp, the other isotopic ratios 
depend dominantly on the circulation rate. The relationship is shown between models 
of CBP as parameterized by a diffusion formalism within the stellar evolution model 
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and those using the mass-flow formalism employed here. They are shown to be effec¬ 
tively equivalent. In general, the CBP treatment readily permits calculation of envelope 
compositions as affected by different degrees of extra mixing, based on stellar structures 
computed by normal stellar evolution models. Using these results, the isotopic ratios 
under conditions of C/0< 1 and C/0> 1 are compared with the data on circumstellar 
dust grains. It is found that the and ^®A1/^^A1 observed for ox¬ 

ide grains formed at C/0< 1 are reasonably well-understood. However, the 
12 c/13c, a,nd ^®A1/^'^A1 in carbide grains (C/0> 1) require that many of their stellar 
sources must have had at least a factor of 4 lower than the solar value. This 

allows a self-consistent description of all these isotopes in most SiC grains. The rare 
grains with 10 cannot be produced by any red-giant or AGB source, nor are 

they reconcileable with novae sources. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: AGB and 
post-AGB — solar system: formation — meteors, meteoroids — stars: abundances — 
Galaxy: abundances 


1. Introduction 

The purpose of this study is to explore the effects of an extra circulation process on the 
^®A1/^'^A1, and ratios in the envelope of an AGB star of solar metallicity. 

By an extra circulation process we mean the transport of matter from the fully convective envelope 
into the underlying radiative region, down to the outer parts of the H-burning shell. The actual 
physical-dynamical basis for such penetration is not known. However, the requirement that some 
such mechanism must sometime occur has been amply demonstrated by observations on RGB stars 
and in studies of dust grains which formed in circumstellar envelopes (e.g. Huss et al. 1994; 
Boothroyd, Sackmann & Wasserburg 1994; Harris, Lambert & Smith 1985; Kahane et al. 1992). 
The isotopic effects of such extra transport on the abundances depend only on the degree of nuclear 
processing. This is directly related to the nuclear reaction rates, the rate of mass transfer and 
the depth to which the extra mixing mechanism penetrates. A detailed physical fluid-dynamical 
mechanism can not significantly affect the resulting abundances beyond determining the rate and 
depth. It is clear that any extra mixing mechanism that does not fractionate one set of elements 
relative to the other elements involved in the nuclear reactions will not change the resulting products 
that are added to the envelope from what is found here. The issue of how small amounts of extra 
mixing will change the stellar evolution is not known but will be commented on. The results, 
presented here with inter-relationships between the different isotopic ratios, are of use in discussing 
and understanding isotopic results obtained by extensive studies of recovered circumstellar dust 
grains. Insofar as a relationship may be found between the isotopic ratios listed above and the 
ratios of major elements in the AGB envelope (G/0 and N/0), this may lead to understanding the 
chemical-physical origins of the grains and the nature of the stars around which they formed. Such 
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results may also be of value in understanding the compositions of planetary nebulae. 

It has been proposed that many short-lived nuclei present in the early solar system (ESS) were 
the products of an AGB star that injected freshly-synthesized nuclei into the interstellar medium 
(ISM) from which the solar system formed (Wasserburg et al. 1994). Because AGB stars are the 
site of s-process nucleosynthesis, there is a rather thorough understanding of the yields of a large 
number of nuclides in their convective envelopes (cf. Busso, Gallino & Wasserburg 1999, hereafter 
BGW99). Among these nuclides, ^®A1 plays a critical role in evaluating the hypothesis that an 
AGB star supplied short-lived nuclei to the early solar system. While it is clear that ^®A1 should be 
produced in AGB stars by the reaction ^®Mg(p, 7 )^® Al (Forestini, Arnould & Paulus 1991; Mowlavi 
& Meynet 2000; Lattanzio, Forestini & Gharbonnel 2000), the quantitative ^®A1 yield, taking into 
account the complexities of destruction of this nuclide, remains very uncertain. As shown by 
Wasserburg et al. (1994) and BGW99, ^®A1 is produced efficiently in the hydrogen shell, but most 
(~ 90%) of this is mixed into the He shell at the onset of intermittent He shell burning. The mixed 
material is then exposed to neutrons in the convective (He-burning) region formed during the He- 
burning pulse, where 80% to 100 % of its ^®A 1 content is consumed by the reaction ^®Al(n,p)^®Mg. 
The neutrons that destroy ^®A1 are provided by the process ^^Ne(Q!, n)^^Mg, activated marginally in 
the relatively high temperatures (T > 3 x 10® K) of the thermal pulse. A value of ^®A1/^^A1~ 10“® 
(= number ratio) results in the envelope after mixing and is much less than the amount initially 
produced in the hydrogen shell. This is a factor of a few less than required if a single AGB star 
produced the Al and this material was then diluted to provide the amount found in the early solar 
system (^®A1/^'^A1~ 5 x 10“^; Lee, Papanastassiou & Wasserburg 1977; Wasserburg et al. 1994, 
1995). Although the discrepancy might be explained by a relatively small change in the fraction 
of ^®A 1 destroyed (i.e., a small variation in the ^®Al(n,p)^®Mg rate), the estimated ^®A 1 yield is at 
best in marginal agreement with the yield required to provide the solar inventory. If the true yield 
is lower than ^®A 1 /^'^A 1 ~ 4 X 10 then an AGB source for solar ^®A1 is not reasonable. There is at 
present no way to improve the estimated ^®A1 yield. Mowlavi & Meynet (2000) recently carried out 
an extensive analysis of Al production in AGB stars. They explored the effects of recently-revised 
reaction rates, and computed the yields expected at varying metallicities and initial masses. The 
results of this work agree with the above expectations: normal, low-mass AGB stars were found 
to be probable sources of ^®A1 abundances observed in most (not all) SiC grains, with ^®A1/^^A1 
values in the same range discussed above. Mowlavi & Meynet also found more efficient production 
of ^®A1 in intermediate mass stars (hereafter IMS) experiencing hot bottom burning (HBB). This 
process was studied extensively by Frost et al. (1998), Lattanzio & Forestini (1999), and Lattanzio 
et al. (2000). However, uncertainties in the final fates of these stars (mainly in the number of 
thermal pulses experienced) prevent a quantitative estimate of their possible contributions to the 
ISM and to the Galactic inventory. 

A direct challenge to standard AGB models arises from the discovery of circumstellar oxide 
grains recovered from meteorites, on which rather precise isotopic ratios are measured. Many 
of these grains reflect ^®A1/^^A1 at their formation far above 10“®. Some oxide grains also have 
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18o/16o about 30% less than the solar value and much greater than the solar value. These 

oxygen data were interpreted to reflect stellar first dredge-up (e.g., Huss et al. 1994, Boothroyd, 
Sackmann & Wasserburg 1994). Studies by Nittler et al. (1994,1997) and Choi, Wasserburg & 
Huss (1999) found some oxide grains with very low or essentially no It was shown that such 
low could result from HBB (Boothroyd, Sackmann, & Wasserburg, 1995), but only in 

stars with masses AMq < M < 8Mq. However, the data on many oxide grains were found to lie in 
a region of the and diagram that is not accessible by HBB. Direct observations of 

oxide molecules in the photospheres of some red giants had for some time indicated similar oxygen 
isotopic compositions, albeit with very large uncertainties (Harris, Lambert & Smith 1985; Harris 
et al. 1987). Similar modifications of oxygen isotopes were subsequently confirmed with a higher 
level of confidence in circumstellar envelopes (Kahane et al. 1992). 

Discrepancies between oxide grain compositions and the compositions that result from simple 
dredge-up in low-mass stars were explained by the hypothesis of extra mixing of material through 
the zone of radiative energy transport above the hydrogen shell - referred to as “deep mixing” or 
“cool bottom processing” (CBP). Such a mixing process is not treated quantitatively in any stellar 
model with a stable and totally non-circulating radiative zone, but it cannot be ruled out on physical 
grounds; this extra mixing has been suggested by many authors (see below). Several abundance 
peculiarities besides the oxygen problem are explained by introducing CBP during at least the 
red giant branch (RGB) phase of evolution. Specifically, observed patterns of ^Li production and 
destruction, and of production on the RGB require GBP (Gilroy & Brown 1991; Charbonnel 
1995; Boothroyd & Sackmann 1999; Gratton et al. 2000). CBP also explains the evolution in time 
of the carbon isotope abundances along the RGB (Denissenkov & Weiss, 1996; Cavallo, Sweigart 
& Bell 1998). The actual mechanism that causes this mixing is not known. It is often ascribed to 
the shear instability in differentially rotating stars (Zahn 1992; Maeder & Zahn 1998), but other 
mechanisms, such as meridional circulation and convective overshoot (Sweigart & Mengel 1979; 
Herwig et al. 1997 and references therein), are available and equally likely. The specific mechanism 
will determine the actual depth and rate of mixing, but should not otherwise matter. 

Deep mixing processes were shown to be especially important in low-metallicity stars, where 
deep mixing may involve relatively hot H-burning cycles. GBP was therefore suggested (e.g. by 
Langer, Hoffman & Sneden 1993, Weiss, Denissenkov & Gharbonnel 2000) as an explanation for a 
collection of CNO, Ne-Na and Mg-Al anomalies observed in low-metallicity stars (reviewed in e.g. 
Kraft 1994, Da Costa 1998). These studies showed that when CBP becomes deep enough to involve 
the H shell itself, the stellar evolution is significantly altered by extensive He enrichment of the en¬ 
velope (Sweigart 1997). Denissenkov & Weiss (2001), investigating elemental abundance anomalies 
in the globular clusters ui Cen and M4, found that at least some of the observed variations in [Al/Fe] 
could result from CBP occurring in RGB stars. These workers used a diffusive model of extra mix¬ 
ing with rate and depth of mixing determined to fit observational trends. They demonstrated the 
possibility of ^®A1 production in RGB stars undergoing GBP, at least in low-metallicity stars and 
with particular parameters. Very recently, observations of planetary nebulae were shown to require 
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further addition of CN-cycle products into the envelope during the third dredge-up phase, even in 
stars with masses that are too low to permit HBB (Pequignot et al. 2000). The destruction of ^He 
that should accompany CBP would also help explain the apparently weak evolution of the Galactic 
^He abundance since the big bang (Wasserburg, Boothroyd & Sackmann 1995; Balser et al. 1998; 
Gloeckler & Geiss 1998; Bania, Rood & Balser 2002). 

For all the above reasons, CBP is now generally accepted as a significant process occurring in 
a second mixing episode late in the evolution of red giants (see e.g. Boothroyd & Sackmann 1999; 
Gratton et al. 2000). The main restriction on CBP is that there be no abrupt increase in mean 
molecular weight anywhere between the convective envelope and the hydrogen shell to prevent 
mixing (Charbonnel, Brown & Wallerstein 1998; Boothroyd & Sackmann 1999). This condition is 
known to be satisfied in M < 2Mq stars on the RGB, at luminosities above that of the “bump” in 
the RGB luminosity function. It should also be satisfied throughout the subsequent AGB phase - 
and certainly during the thermally-pulsing AGB (TP-AGB) phases - for all low-mass stars. From 
all these considerations CBP therefore appears to be both possible and necessary in the evolution 
of low-mass stars. 


2. The Model 

The general structure and evolution of AGB stars have been studied extensively (for a thorough 
introduction see Lattanzio 2001^). The structure of a low-mass star during this phase of evolution 
is schematically as follows (see Fig. 1): in the center is a degenerate carbon/oxygen core. Outside 
this is a helium region which cannot support steady nuclear burning. The next layer is a thin 
(< 10 “^Mq) region where hydrogen burning takes place, providing essentially all of the star’s 
luminosity. The outward advance of the H shell adds to the He region, which periodically undergoes 
nuclear burning (“thermal pulses”) in short (~ 200 yrs) bursts separated by long (~ 30,000 yrs) 
intervals. The H shell and the region immediately above it (< 10~^Mq in total) are stable against 
convection, and we call these together the “radiative zone.” Finally, the outer region of the star is 
a large, fully convective envelope. When mass loss is included as described below, this envelope has 
mass ~ 0.9Mq at the start of the TP-AGB and shrinks as the H shell advances outward and as the 
star sheds envelope material in a wind. Late in the AGB evolution, each thermal pulse is followed 
by a brief disruption of the structure just described, during which material that has seen nuclear 
processing in the H and He shells is mixed into the envelope. This mixing is called “third dredge- 
up” (see review by BGW99). In the above description, the radiative zone contains no circulation. 
We will treat the problem of extra mixing following Wasserburg, Boothroyd & Sackmann (1995; 
hereafter WBS). Here a slow circulation is added by imposing transport of envelope material at 
a mass flow rate M down to some depth inside the radiative zone specified by temperature Tp 
and then back to the envelope. Typically, a small amount of nuclear processing takes place in the 


^John Lattanzio’s online stellar evolution tutorial available at http://www.maths.monash.edu.au/ johnl/stellarevoln/. 
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circulating material. 

We use a numerical model of an AGB star with solar initial composition and I.SMq initial 
mass. The model was computed by Straniero et al. (1997) using the FRANEC code to follow 
the full evolution of the star from the zero-age main sequence through the TP-AGB phase, with 
no mass loss. The corresponding loss of envelope mass was subsequently computed using the 
Reimers (1975) prescription as discussed by Busso et al. (2001). The choice 77 = 1 for the free 
parameter of the Reimers formula was adopted. Details of the mass loss and dredge-up are relatively 
unimportant for our conclusions: runs of the code at a fixed envelope mass of 0.5Mq produced the 
same results to ~ 10%, and variations in dredge-up efficiency will be discussed below. The detailed 
evolution on the AGB of the above model is summarized in Fig. 2 and Table 1. We will refer 
to the detailed structure and evolution of this model, including both dredge-up and mass loss, as 
the stellar evolution model, SEM(1.5). For our calculations, we assume that extra mixing may be 
modelled by studying circulation inside the radiative-zone structures previously computed in the 
SEM(1.5) model. The nuclear processing in GBP depends mainly on the variation of temperature 
T and density p as functions of position within these radiative-zone structures, and of time t. These 
structures are the environment within which the circulating material moves and is processed. The 
SEM(1.5) model also provides us with the timescales for AGB evolution, the rate of mass loss 
from the envelope, and the changes to envelope composition due to third dredge-up (TDU). The 
SEM(1.5) model exhibits TDU after the 10th thermal pulse, and after each subsequent thermal 
pulse until the 26th (the interval from Stage 5 to Stage 6 in Fig. 2 and Table 1; two subsequent 
pulses produce negligible dredge-up). This dredge-up brings large quantities of newly-produced 
carbon and of s-process nuclides to the envelope, and as a result, the envelope of SEM(1.5) with 
no GBP has C/0> 1 from the 22nd pulse onward. The envelope begins GBP with a composition 
reflecting earlier first dredge-up, as discussed below for individual nuclides. 

It is useful in discussing GBP to distinguish between three types of material. “Envelope 
material” (E) makes up the convective envelope and has abundance (in moles per gram) of 
nuclide i. “Circulating material” (C) is that material moving through the radiative zone. It starts 
out as a representative sample of the envelope at a given time, and then moves down through the 
radiative zone and back as shown in Fig. I. Abundances Yf' in the circulating material are changed 
by {p, 7 ) and {p, a) reactions, as well as beta decays along the flow path. The circulating material 
just before returning to the envelope is the “processed material” (P), with composition denoted 
Yj^. We assume that mixing of processed material with the convective envelope is instantaneous. 

The numerical GBP model follows the envelope abundances as material descends into the 
radiative zone and is replaced by newly-processed material. This corresponds to integrating 


Y^ = 


M 

Me 




( 1 ) 


as, Me, Y/’, and evolve together. The quantity Y/^ is computed by following the position of 
a parcel of material as it starts at the bottom of the envelope, circulates downward through the 
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SEM(1.5) stellar structure, and finally returns to the envelope with composition Yf. The stellar 
structure and M specify as functions of time the conditions of T and p seen by this material, and 
the numerical code integrates a nuclear reaction network (rates from Angulo et al. 1999; discussed 
in Appendix A) with these changing conditions. We do not calculate the detailed evolution by 
following every detail of development of the radiative zone with time. For a given calculation, we 
choose one of the times from the SEM(1.5) evolution listed in Table 1 (we call this a “stage”), 
and we hold the radiative-zone structure constant as specified by this stage. This amounts to an 
approximation that the processing is independent of the evolution of the radiative zone for the 
duration of CBP. A check of this assumption is given in Sec. 3. We consider the dominant extra 
processing to take place when the H shell temperature Th (defined as the temperature where the 
maximum energy generation takes place) exceeds logTn > 7.7 and when TDU is active. This is 
restricted to the last 2 x 10® yrs of AGB evolution of the I.SMq star. When considering CBP during 
the interval between Stages 5 and 6 , we use the radiative-zone structure from Stage 5, and we include 
TDU by interrupting the evolution at the times of the thermal pulses and altering the envelope 
abundances to reflect dredge-up as computed for SEM(1.5). The CBP code is re-started after each 
pulse with the new composition as the intial state. We refer to the results of this calculation as the 
“numerical model.” We will also present analytical arguments to interpret the output. After an 
initial examination of the conditions necessary for ^®A 1 production and incorporation into grains, 
we confine ourselves to evolution from Stage 5 to Stage 6 , when TDU is active, and we use the 
radiative-zone structure from Stage 5. All of the results from the numerical model shown in the 
figures correspond to these calculations. 

A fully self-consistent model of CBP using this “conveyor belt” model would involve calculating 
the abundances of the nuclides in the envelope at all times and following the evolution of stellar 
structure as conditions in the radiative zone change, as the H-shell advances into regions that have 
been processed by CBP, and as dredge-up and mass loss change the envelope. The most extreme 
assumption that we have made is that a radiative-zone structure from a fixed time provides a 
reasonable approximation to the radiative-zone structure for the next few 10® yrs. We find below 
that this is true during the TP-ACB (excluding very brief interruptions by thermal pulses), and 
that C, N and O isotopes are insensitive to the specific radiative-zone structure chosen. We will 
show that for ^®A 1 there is a strong dependence on this structure, dominantly due to the evolving 
Th- It will be shown that there are nonetheless clear general rules for ^®A1 that adequately describe 
the system. Note that the conditions of T and p in the envelope play no role in this problem as no 
nuclear processing occurs there. 

The radiative zone is only 8.6 x thick at Stage 5, so if the envelope has a convective 

overturn time of ~ 1 yr, and the radiative zone turns over much less frequently than the envelope, 
then M < lO“^M 0 /yr, roughly. The other important constraint is that stellar energy generation 
not be disrupted too much, or approximately that Tp < Tp- In principle, M does not completely 
specify the speed of matter circulation. The speeds of the upward and downward streams may be 
expressed as in WBS by fu and fd, respectively the fraction of the area at fixed radius occupied by 



the upward and downward streams. We have fixed both parameters at fu = fd = f = 0.5, and they 
will not have any significant effect for fixed /„ + fd (see WBS). For each calculation we fix a value 
of M and a maximum depth of penetration in the radiative zone corresponding to a temperature 
Tp. There are no other independent parameters in the calculation. 


2.1. 26^1 

We first consider ^®A1 which is the simplest case. The long-lived ground state of ^®A1 is 
produced in the circulating material by the process ^^Mg(p, 7 )^®A 1 and destroyed by ^®Al(p, 7 )^'^Si. 
The ^®A1 abundance in a sample of circulating material starts with the current composition 
of the envelope and then evolves during the time that it circulates through the radiative zone 
according to 

5^26 = YgX[^^Mg{p,j)] - YgX[^^Al{p,^)], ( 2 ) 

where A[rxn] = pNA{c^v)Yg for the reaction labeled “rxn”, p is mass density, Na Avogadro’s 
number, {av) the thermally-averaged reaction rate, Yg the ^^Mg abundance, and Yg the hydrogen 
abundance. The quantities on the right-hand side of this expression evolve as the material’s position 
in the radiative zone changes. The destruction term is negligible for two reasons; 1) for the 
presently-recommended rates (Appendix A and Angulo et ah 1999), the destruction cross section 
is two orders of magnitude smaller than the production cross section, and 2) Yg never builds up 
to a level comparable to Yg. The rate of ^®A1 decay during the time spent circulating through the 
radiative zone is negligible. 

For most values of M and for Tp < Tp, the ^®A1 production rate is low enough that only a 
small fraction of the ^^Mg burns on a single pass through the radiative zone. As a result, Yg is 
almost constant, and the amount of ^®A1 in the processed material is approximately 

>^26 (i/) = Yi{U)+ f'YgX[^^Mg{p,j)]dt (3) 

Jti 

fMgCE 

y 26 itf) = Yi{ti) + 2fM-^YiYg pNA{av)dm, (4) 

JMp 

where we have taken advantage of the relationship between M and the time-dependent position of 
the circulating matter to write the integral in mass coordinates with no explicit reference to time. 
Mp is where T = Tp and MpcE is the bottom of the convective envelope (see Fig. 1) - the time 
integral begins at the time ti when the circulating material leaves the envelope and ends at time 
tf when it returns to the envelope. (The factor of 2 arises because the circulating material passes 
each point twice: once on the way down, and once on the way back up.) Note that Yg{tf) — Yg{ti) 
does not change with time, except to the extent that stellar evolution changes the path integral. 

The abundance of ^®A1 in the envelope reflects the addition of processed material with the 
^®A1 abundance Yg to the envelope of mass Mp, competing against free decay with mean lifetime 
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T26 = 1.0 X 10® yr and removal of material into the downward stream: 

'E _ M , p 

-'26 — I -* 26 -'261 "26 -'26 

i-Mbce 

= 2/M/yJ(0)y# / pNA{av)dm-T^^^Yi 

JMp 

— “ ^26 -' 26 ) 


(5) 

( 6 ) 
(7) 


where the last two lines define the production rate P per unit mass of the envelope. Note that 
all factors of M have cancelled: the amount of processed material is proportional to M, but the 
time it spends at high temperature is proportional to M~^. If the envelope mass decreases with 
time and the reaction path integral over the circulation stays constant, then P will increase with 
l/ME{t). 

Calculations were done in the numerical model using the radiative-zone structure of Stage 5, 
with CBP occurring roughly for the duration of third dredge-up (from Stage 5 until Stage 6). The 
results for ^®A1 are shown in Fig. 3a as a function of Tp. During this interval, ^®A1/^^A1 = 6 x 10“^ 
is added gradually by normal third dredge-up, and this amount is shown as the “TDU only” line. 
As argued above (cf. Eqn. 6), there is no significant dependence on M. We see that for log Tp < 7.6 
(0.2 dex below Tp) there is no significant ^®A1 production. The term log (^®Al/^^Al) rises almost 
linearly with log Tp as Tp increases toward Tp- It is evident that CBP penetrating to temperatures 
above Tp = 10^ ®K during this phase of evolution can readily produce ^®A1/^^A1 as high as found 
in extreme examples of circumstellar oxide grains. 


In order to understand these numerical results, we consider analytic solutions to the approx¬ 
imate expression in Eqn. 6. Table 1 and Fig. 2 show the selected stages of AGB evolution on 
which we will now focus. Time zero is taken as the state when He core burning is just finished, 
and qualitative descriptions of the evolutionary stages are given in Table 1 along with character¬ 
istic parameters of the radiative zone. For the semi-analytic treatment, we neglect mass loss and 
dredge-up, and we assume that the stellar structure (i.e., the structure at a selected stage from 
Table 1) does not change with time. Using the solution to Eqn. 7 for a constant-mass envelope 
and a selected time-independent radiative-zone structure, we obtain 


Yi = yj (0)e-*/”^s + 


1 — e 


-t/T26 


T2&P- 


( 8 ) 


Let us now consider the steady-state abundance of A1 that could be achieved in an AGB star 
fixed at the structure from a particular evolutionary stage. The results for Stage 5 are shown in 
Eig. 3a and are in good agreement with the corresponding numerical calculation. We consider the 
following additional times in the evolution: 1) The stellar structure at the immediate end of core 
He burning; 2) the middle of the early AGB; 3) just before the first thermal pulse; 4) just after the 
second mild pulse, etc.; and finally 6) just before the twenty-sixth thermal pulse (see Table 1 and 
Eig. 2). The ^®A1/^’^A1 ratios at steady state {t 2 &PjY,^) for each of these stellar stages are shown 
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in Fig 3b. There is no dependence on M because of the cancellation noted above. The end points 
of the curves correspond to the extreme limiting case of Tp = Tp. By the onset of the first thermal 
pulse, the production rates rise rapidly and ^®A1 production is possible only at these late stages. 
The critical parameter is Tp, which is not a priori known. Values of log Tp > 7.6 are required to 
produce substantial amounts of Al. Inspection of Fig. 3b shows that to achieve the same value of 
^®A1/^'^A1 requires different in Tp for each choice of stellar structure (e.g., ^®A1/^'^A1= 0.03 requires 
logTp = 7.64 at Stage 3, 7.68 at Stage 4, 7.72 at Stage 5, and 7.73 at Stage 6). However, if we had 
chosen Tp/Tp ~ constant for each of the last three stages rather than fixing Tp, then the ^®A1/^'^A1 
production rate would be approximately constant. (See section 3 below.) 

In general, high levels of ^®A1/^^A1 can be achieved by CBP, essentially independent of M, 
but completely dependent on Tp. Thus, ^®A1/^'^A1 is a measure of the depth of penetration, 
independent of M. Significant ^®A1 production can occur over all the late stages, but we expect it 
to be dominated by production over the last T 2 q of time before the envelope is lost. It can be seen 
in Fig. 3a that the numerical model is very close to the simple analytical calculation. Most of the 
small discrepancy arises because the TDU period lasts only for 8.5 x 10^ yr, not long enough to 
reach steady state. Overall, the agreement is astoundingly good. 

In summary, significant ^®A1/^^A1 can result from CBP late on the AGB. The exact amount 
of ^®A1 produced depends on the depth of circulation as expressed by Tp and the duration of 
CBP (limited by the AGB lifetime), but only very weakly on the circulation rate M. To produce 
26A1/27A1 > 4x 10 ^ just before the star ejects its envelope, it appears necessary that logTp > 7.67. 
There is only small ^®A1 production in the early AGB phase and we infer that this also applies to 
any CBP during the preceding RGB phase because of the low Tp at that time. 


Table I: Stages of AGB evolution considered for CBP. “Start time” is in years after the end of core 
He burning. AM is the thickness of the radiative zone starting at the bottom of the H shell. The 
numerical model is based on Stage 5, assuming a constant structure for the radiative zone for the 
interval until Stage 6. 


Stage 

Start Time (Myr) 

Description 

logTp 

AM {Mq) 

Me (Me) 

I 

0 

end core He burn 

7.423 

0.376 

0.587 

2 

21.04 

before H-shell re-ignition 

7.456 

0.0329 

0.900 

3 

22.47 

before 1st pulse 

7.679 

5.18 X 10-3 

0.862 

4 

22.78 

after 2 mild pulses 

7.740 

2.02 X 10-3 

0.887 

5 

23.66 

before lOth pulse (start TDU) 

7.797 

8.65 X 10-^ 

0.741 

6 

24.49 

after 26th pulse (end TDU) 

7.834 

3.17 X 10-^ 

0.242 
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2.2. Oxygen 

We now examine the implications for the oxygen isotopes of the conditions that produce ^®A1. 
In contrast to ^®A1, is not produced but is destroyed. First dredge-up during the RGB phase 
provides the envelope with that is decreased from its initial value due to destruction during 

the main-sequence phase. This has been shown by Dearborn (1992), who recognized the significance 
of and as indices of interior stellar temperatures and mixing processes. These results have 
been confirmed by El Eid (1994), Boothroyd & Sackmann (1999), Gallino et al. (1998) and others. 
In contrast, the isotope is produced in the interior and is then enriched in the envelope above 
the initial value as a result of first dredge-up. Neither nor is affected by thermal pulses. 
The major isotope in the envelope, is essentially unchanged throughout stellar evolution. 

The results of the numerical model for and are shown in Eigs. 4 and 5 as 

functions of M. The curves show the final envelope composition after using the fixed radiative-zone 
structure of Stage 5 and processing through the length of the TDU period. Curves for various values 
of logTp are given. It can be seen that for M < lO“^M 0 /yr, there is no substantial change in 
for all Tp. This is because very little material is processed over the lifetime of the AGB for these 
mixing rates. Eor M > 3 x lO“^M 0 /yr there is a sharp drop-off with a magnitude depending on 
Tp. For rates above 3 x lO“®M 0 /yr, the destruction is almost complete when log Tp > 7.62. The 
horizontal parts of the curves at higher Tp are due to a small amount of production balancing 
the destruction. For there is a rise to higher by a factor of ~ 1.5 (dependent on Tp) 

above the value after first dredge-up. We find that if log Tp exceeds a low threshold, the calculation 
of the final composition as a function of M does not change significantly. This threshold is 7.62 at 
Stage 5 and 7.60 at Stage 3. Therefore, the results are robust and a calculation of the full evolution 
with continuously changing structure and GBP will not significantly alter the results found here by 
assuming a fixed radiative zone. 

In order to understand the results, we examine the transport equations as was done for ^®A1. 
The basic results can be obtained almost quantitatively from an analytic treatment. As long as 
logTp > 7.62, the evolution of in GBP is characterized by rapid destruction. (This applies for 
Stage 5; other stages are characterized by analogous thresholds, never lower than 7.60 or higher than 
7.65.) This destruction occurs via ^®0(p, a)^^N, with abundance in the circulating material 
evolving according to 

yH = -X[^^0{p,a)]Yg. (9) 

(Production via ^^0(p, 7 )^®F(/ 3 '’'p )^®0 is negligible). The abundance in the processed material 
returning to the envelope is the time integral of this equation along the entire path taken by the 
material through the radiative zone, given by 

^18 = A[^®0(p,a)]dt| (10) 

f pMeCE 

= Eifexp - 2 /M-iy^ / 

I JMp 


pNA{crv) dm 


( 11 ) 
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The exponential factor is typically very small for even modest values of Tp; for example, this 
expression gives = 0.007 if log Tp = 7.67 and M = 10“® Mq/ji using the radiative zone 
of Stage 5. The result is as low as 10“^ in the processed material. Higher Tp or lower 

M would further reduce as would processing slightly earlier in the TP-AGB evolution when 
the threshold for destruction is lower (analogously with ^®A1 production). For example, using the 
radiative zone of Stage 5 with logTp = 7.69 and M = 10“® Mq/yt gives Y^^ = 4.2 x 10“^Y^f. 
Using the radiative zone at Stage 3 with M = 10“® Mq/jt, only logTp = 7.61 is needed to produce 
Yif = 0.007Yif. 

Because destruction in the processed material is so thorough, the envelope abundance of 

reflects only the removal of envelope material and its replacement with material containing 
essentially no 


- -Me 


The result is exponential decay of with time, so if Me is constant, 


n 


%pexp(-Mt/Mp). 


( 12 ) 


(13) 


In this analytical treatment, the destruction of in the processed material is taken to be 
almost complete. This will apply to all stages of the evolution above the threshold Tp discussed 
above. Hence the resulting in the envelope at the end of the AGB is given by Eqn. 13 

where the only parameters are M,Mp, and the time scale for AGB evolution. This is shown for 
t = 8.5 X 10^ yr (the length of the TDU period) and Mp = 0.5 Mq as the dashed curve in Figure 
4. It can be seen that the semi-analytic treatment is almost indistinguishable from the numerical 
model over most of the range of interest. The flat region at M > lO“^M 0 /yr would have been 
obtained if we had not neglected the small production. 

Because of the low temperature at which is destroyed, it is evident that low in 

the envelope is not necessarily related to ^®A1 production. As long as Tp > 7.62, in the 

envelope is completely independent of Tp. This is illustrated by plotting in Fig. 6 for 
vs. ^®A1/^'^A1 the final compositions after GBP running through the TDU period. Compositions 
were computed with the numerical model. One can see that the curves of constant M (dotted 
lines) are approximately curves of constant final while the curves of constant Tp (dashed 

lines) are approximately curves of constant ^®A1/^^A1. The departure of the constant-M curves 
from constant is explained as follows: at low Tp (low ^®A1/^^A1), destruction on a 

single pass through the radiative zone is incomplete; at high Tp (high ^®A1/^^A1), the situation is 
complicated by partial destruction. 

The final envelope compositions are the compositions most likely to be reflected in presolar 
grains; the most grain formation should occur with very high mass loss at the very end of AGB 
evolution (Table 1). The circles in Fig. 6 show the abundances found in oxide grains. It appears 
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that the model can, in principle, explain simultaneously the ^®A1 and contents of all the oxide 
grains found so far, and with reasonable values of M and Tp (however, see the trends below). 
The largest value of ^®A1/^^A1 indicates a maximum Tp in the source stars of about K, while 

a wide range in M is suggested if all the values are attributable to CBP. Not all grains 

formed from 1.5 Mq stars. Moreover, inferred values of M and Tp are sensitive to exactly when 
CBP occurs and for how long; these values do provide qualitative insight into the actual mixing 
rate and depth. In any case, decreases at a rate dependent upon M and not Tp, so that 

there is no strong correlation with ^®A1/^^A1 production except that ^®A1/^^A1 and may 

be modified simultaneously by CBP. Low M and high Tp will produce small depletions and 
abundant ^®A1. High M and high Tp will produce major depletion and high ^®A1. Thus, 
a direct correlation between ^^Al/^'^Al and should not be found, although Nittler et al. 

(1997) and Choi, Wasserburg & Huss (1999) have suggested that it might occur on an empirical 
basis. 

The plateau of at M > 10“®MQ/yr for all Tp seen in Fig. 5 shows net increases above 

the dredge-up value ranging from 0% to 40%. This plateau value increases with temperature, but is 
at 0.0011 for all Tp > 7.62. Fig. 7a shows vs. for the final envelope composition, 

at various values of M and Tp. Generally speaking, the results lie in the space between point D 
(composition after first dredge-up) and the end point B. If we use Stage 3 instead of Stage 5 to 
specify the stellar structure, this changes the results by shifting the end point B to 0.0016. 

Results for Stages 4 and 6 are similar to those of Stage 5. 

We now examine using the semi-analytic treatment. In the circulating material, 

destruction is very small, and the abundance is governed by competing production via 
^®0(p,7)^’^F(/3 +p)^'^ 0 and destruction via ^’^0(p, a)^"^N: 

yfj = 1^i?A[160(p,7)] - YgX[^^0{p,a)]. (14) 


Starting from a low initial abundance, increases as matter moves through the radiative zone 
until the two terms on the right-hand side of Eqn. 14 cancel. The equilibrium ratio specified by 
this cancellation is determined completely by reaction rates so that processed material has 


yP 

ill 

YP 


A[^^0(p,7)] 

A[i70(p,a)] 


~ 0 . 0011 . 


(15) 


As a result, 

^17 ^17 -MtjME 

uf ~ u? (uf(o)jj 


( 16 ) 


The final envelope composition using this approximation is shown as a function of M in the dashed 
curve in Fig. 5. As with the numerical results presented, this assumes the radiative zone of Stage 
5 and t = 8.5 x 10^ yr. We chose Mp = 0.5Mq for the analytic approximation. It is clear that 
at the condition of logTp > 7.62 required to obtain equilibrium, the results are in close agreement 
with the numerical calculation. Thus, the envelope as a function of time reflects the 

replacement of envelope matter with material that has 0.0011. 
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The steady-state is almost temperature-independent over the temperatures charac¬ 

teristic of the lower part of the radiative zone late on the AGB, as shown in Fig. 8 . As a result, 
Y^/Y^ is approximately 0.0011 for all cases with log Tp > 7.62, with an uncertainty of about 30% 
in the ratio arising from measurement errors on the cross section inputs (Blackmon et al. 1995; 
Angulo et al. 1999; see also Appendix A). We note that only in the last six years have 
cross sections been sufficiently well-measured to allow us to reach this conclusion, because of the 
measurement of the 66 keV resonance by Blackmon et al. (1995). At lower temperatures, the 
equilibrium increases, but the timescale to reach equilibrium also increases so that these 

higher values of are inaccessible after Stage 4. The highest value that we found at any 

time on the TP-AGB is 0.0016. 

Since the envelope as a function of time reflects gradual replacement of envelope 

material with material that has fixed = 0 . 0011 , and the processed material also has 

extremely low we now recognize the numerical results in Fig. 7a as mostly falling on a 

two-component mixing curve. 

We also carried out calculations using radiative-zone structures from other stages, and the 
results were almost indistiguishable as long as log Tp > 7.60. As explained in the analytical 
treatment, the dominant effect for logTp > 7.62 is to produce an envelope that is a mixture of the 
envelope composition after first dredge-up and the end point B that represents the equilibrium value. 
For all values of logTp > 7.62, the numerical results lie on the mixing line at positions depending 
almost exclusively on M as indicated in the figure. For lower temperatures, (logTp < 7.62) the 
final compositions fall in a space bounded on the one side by the mixing curve and on the other by 
the curve with M = lO“^M 0 /yr. More generally, as the oxygen isotopic evolution resulting from 
GBP depends on the value of and after first dredge-up and not on the detailed 

stellar evolution, we may readily infer the results of GBP for low-mass AGB stars of generally 
solar chemical abundances but with arbitrary oxygen isotopic composition: the end point B is the 
natural final composition for material that has been near the H shell in any AGB star of low mass 
(M < 3Mq), because all of these stars have H shell temperatures near 10^'® K. 

The results for two initial compositions are shown in Fig. 7b. The upper point A' represents 
ratios in the envelope at the end of first dredge-up for some star. For logTp > 7.62, the trajectory 
is a mixing line between A' and B, where the position only depends on Mt. If logTp is substantially 
less than 7.62, but greater than 7.53, it will be a curve that drops from A' (due to destruction) 
and then swings toward B as M increases. This same rule applies to A!' or any other point, insofar 
as the initial stellar mass does not exceed 4 - 5 M 0 , where HBB may occur (Iben & Renzini 1983). 
As the initial oxygen isotopic composition of the star which produced a particular oxide grain is 
not in general known, we will use this construction in discussing the data on grains in Section 4.1.1. 

In light of the above considerations, the production of 0.0011 by first dredge-up is 

explained because material in the stellar mantle is processed during the main-sequence phase at 
much lower temperatures than those present on the AGB. The relevant reaction rates are much 
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slower and the equilibrium ratio larger (see Fig. 8 ), while the duration of the main-sequence phase 
is much longer than that of the AGB phase. 


2.3. Carbon 

The evolution of carbon in CBP is more complex than that of ^®A1, or since the 
nuclides discussed so far are negligibly affected by TDU. The abundance of in the envelope 
is altered at each dredge-up episode. Consequently, the ratio is shifted. The change in 

C/0 in the envelope at each thermal pulse in SEM(1.5) is shown in Figure 9a along with the 
change in the mass of the envelope with time. In the numerical model, we assume that CBP starts 
immediately after the first thermal pulse that produces dredge-up (Table 1 and Fig. 2). This is 
the same calculation used for ^®A 1 , and 

Figure 9a shows the evolution of C/0 for SEM(1.5) as a function of time up until the end of 
TDU and includes the cases of CBP with M = Mq/yt and log Tp = 7.696 for comparison. 
Note that C/0 = (^^C -|- ^^C)/^®0. In Eigure 9b we show the evolution of ^^C/^^C in the envelope 
versus ^®A1/^'^A1 with time for three cases. As already stated, it is evident that high values of 
Tp give high ^®A1/^^A1. However, to obtain low ^^C/^^C at high ^®A1/^^A1 requires high M (e.g. 
IO“^M 0 /yr). In the latter case, the ^^C content of the envelope is so low that every thermal pulse 
gives a spike in ^^C/^^C which is then destroyed by rapid CBP. The final envelope compositions in 
terms of C/0 and ^®A1/^'^A1 are shown in Figure 10. It can be seen that the approximate boundary 
for C/0 > I is M ~ I0“®'^MQ/yr. The relationship between ^^C/^^C and C/0 is shown in Figure 
II. The solid triangle represents the final envelope composition of the AGB star computed by 
SEM(I.5) with no CBP. The results for CBP are shown as curves of constant M with the values of 
logTp varying with position along these curves. It can be seen that to obtain both C/0 > I and 
low ^^C/^^C ( < 20) requires conditions that are very restrictive ( log Tp < 7.616). Low ^^C/^^C 
may readily occur for higher Tp but with C/0 < I. 

In the material circulating through the radiative zone, the ^^C and ^^C abundances are gov¬ 
erned by the equations 

^2 = Ui^A[i5N(p,a)]-yi^A[i2c(p,7)] (17) 

Yg = YgX[^^C{p,Y)]-YgX[^^C{p,Y)]. (18) 

As a sample of circulating material descends into the radiative zone and its temperature rises, the 
negative term on the right-hand side of Eqn. 17 grows and ^^C is destroyed. The ^^C that is 
destroyed is converted to ^^C, so that quickly both the production and destruction terms in Eqn. 
18 balance. The condition that ^^C/^^C in this material does not change with time then gives the 
steady-state ratio 

Yg _ A[i3C(p,7)]-Apc(^^^)] 

Yg Apc(p,7)] 


(19) 
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which depends only on cross sections. For the temperatures in the radiative zone, it is about 3. 
As processing continues, continues to burn while remains approximately constant, 

producing Finally, if Tp is sufficiently high, enough is made to close the CN cycle with the 
reaction chain ^^N(p, 7 )^® 0 (/ 3 '*'z/)^^N(p, q;)^^C, and the amounts of and become steady with 
I2c/i3cr^j 4 and ~ 0.03[Yj^f(t) + Y^{t)]. Thus, CBP destroys mixing ^^C-depleted and 
^^C-enhanced material into the envelope. This explains the strong dependence of the C isotopes on 
M: because at most temperatures, M sets the rate at which processed material of the new 
composition is mixed into the envelope. 

The effects of CBP can be seen as follows for a given M. The ratio first decreases 

with increasing Tp due to increasing input of processed material, approaching 3. For 

log Tp > 7.6, most of the C is converted to in the processed material, so that processed material 
contains only small amounts of and As a consequence, in the envelope remains 

near the value it would have without CBP. For M > lO“®'^M 0 /yr, the processing at log Tp > 7.6 
is so effective that a substantial amount of carbon is converted to N and C/0 becomes less than 
one. 


2.4. C-rich vs. O-rich envelopes 


The conversion of a large fraction of the carbon in the stellar envelope into nitrogen has impor¬ 
tant observational effects: the C/0 ratio determines whether the signs of related nucleosynthesis 
should be expected in S stars versus carbon stars, or in oxide grains versus carbon and carbide 
grains. As seen in Figs. 9a and 10, when (and whether) this ratio exceeds unity depends on com¬ 
petition between third dredge-up and CBP in the invervals between dredge-up events. To look at 
this competition, we now examine the problem analytically. 

The rate of addition of ^^C to the envelope by TDU is some function of time, F{t). CBP 
removes material from the envelope and replaces it with material in which some of the ^^C has 
been destroyed. The time evolution of the envelope ^^C abundance is then governed by the equation 


^ (^12 - ^12) + (20) 

We wish to find the condition on M such that the envelope is carbon-rich at the end of the TP-AGB. 
IIYP = {1 — a)Fi 2 (where a is the fraction of ^^C destroyed in a single pass), we then have 


Yi = -—aYi + F/ME. ( 21 ) 

F{t) is essentially a sum of delta functions from dredge-up, but we will approximate it by its mean 
value, {F) = 12 x 4.3 x lO^^T^f (0)M£; yr ^ (see Fig. 9a). The solution to Eqn. 21 is 

Tjf (t) = (0) {1 - . (22) 

Ma ^ ) 
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At the end of the TDU phase, t = 8.5 x 10^ yr. The condition that the envelope is carbon-rich is 
that Y^/Y^ > 1. Inserting this condition along with Y^{Q) = 1.7 x 10“^ into Eqn. 22, dividing 
by Y^ = 5.9 X 10“^, and solving for M gives 

M < 7 X Me X 10“'^ yr“^/a. (23) 

The value of the fraction destroyed (a) depends strongly on Tp and weakly on M. A reasonable 
mean value of Mp during the TDU period is 0.4Mq. If a 1, then M < 3 X 10 for 

C/0 > 1, within a factor of a few of what we find numerically. If M = lO“^M 0 /yr, the maximum 
CBP rate we considered, then C/0> 1 in the envelope requires a < 3 x 10“^, or very little carbon 
destruction per pass. The combination of high M and C/0> 1 is therefore restricted to very low 
Tp. 

If the true efficiency of dredge-up is different from that for SEM(1.5), then (F) differs from 
the value used here (for a possible mechanism, see Herwig et al. 1997). If (F) is doubled, then 
the M that divides C/0 > 1 from C/0< 1 (Eqn. 23) increases by a factor of three; we find that 
the relationship between the two parameters is approximately linear over this range of {F). Since 
dredge-up does not have a strong effect on the other abundances, the only influence on our results 
is to change which M lie on which side of C/0< 1. The abundance changes found above for the 
O and below for the N isotopes take place over wider ranges in M than factors of three, so we do 
not expect qualitative changes in our results if the dredged-up compositions change. 


2.5. Nitrogen 

The effect of CBP on nitrogen is to efficiently destroy and to greatly increase in the 
envelope. The increase in is due to the conversion of ^^C into by CBP at log Tp > 7.55. 
Complete conversion of ^^C in the envelope from the initial C inventory in the absence of TDU 
would increase the abundance to over three times the original main sequence value. As TDU 
considerably increases the amount of available carbon in the envelope, even more will be 
produced. This can, at the end of the stellar lifetime, give 1.5 < N/0 < 8.6 (with values above 
1.7 involving destruction of oxygen). Such large amounts could constitute an important source of 
primary (scaling linearly with metallicity) nitrogen in the Galaxy. 

The relationship between and ^®A1/^^A1 is shown in Pig. 12a for the final envelope 

composition as a function of Tp and M for the numerical model. It is seen that there is a rapid 

drop in with increasing M for all logTp > 7.50. For logTp > 7.62, the ratio is 

essentially constant for a given M. Comparing with Fig. 10, we see that a wide range in 
(from about 3 x 10“® to 6 x 10“^) may occur if C/0> 1. If C/0< 1 then M > 10 “®'^MQ/yr and 

and 3 x 10“^. See the discussion of SiC grains in Section 4.1.2. 

In Fig. 12b we show N/0 vs. C/0 for the final composition of the envelope in the numerical 
model. The curves for different M and Tp are indicated. It can be seen that CBP will produce high 
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N/0 for M > 10 and log Tp > 7.6. We note that C/0> 1 requires M < 10 ®'^MQ/yr 

for the same Tp. 


3. Self-consistency and relation to other models 
3.1. Effects on energy generation 

We now examine the self-consistency of the post-processing approach with the basic AGB 
stellar model. Since we assume an underlying stellar model [SEM(1.5)] without CBP, it is important 
that the additional processing not affect the energy generation in a way that would significantly 
alter the evolution. The important issue is to determine what values of Tp and M plausibly allow 
a self-consistent calculation. 

For M, we have considered the range 10“® to lO“^M 0 /yr. The low end is comparable to the 
rate of advance of the H shell, and too slow to turn over the radiative region between two successive 
thermal pulses. At the high end, the rate at which the radiative zone turns over corresponds to 
M/AM ~ 0.1 times per year. If the CBP circulation is to be distinct from envelope convection, 
its turnover rate should be at least an order of magnitude less than that of the convective envelope 
(~ 1 per year). Concerning the effective penetration depth Tp, we chose the lower limit of 7.44 for 
negligible processing and the upper limit of 7.76 to produce ^®A1/^^A1~ 10“^. With this range in 
parameters we now examine the effects on the energy production rate and luminosity. The energy 
generation in the parts of the stellar interior where CBP takes place occurs at very nearly the same 
rate with and without CBP. This is because CBP never changes the H mass fraction on a single 
pass through the radiative zone by more than 10“^ in the CBP models. It was calculated that 
CBP will change the H mass fraction in the envelope over the whole evolution by a total of only 
4.6% even at the highest Tp we considered. Since we have assumed that CBP changes neither 
the physical conditions in the radiative zone nor the total number of C and N catalyst nuclei, the 
energy generation rate can change by at most 4.6%. An estimate of the energy generated in the 
matter undergoing CBP is computed as follows; we examine the H mass fraction of the envelope 
at the end of the numerical calculation and find the difference from the final mass fraction in the 
original SEM(1.5) model in the absence of CBP. We multiply this difference by an average envelope 
mass of 0.6Mq and then by 7 MeV/nucleon, and divide by the 8.5 x 10^ yrs that CBP was active. 
This yields an effective luminosity, shown in Fig. 13 as a function of Tp. For M > 10“^MQ/yr, 
this luminosity is independent of M (because, just as in the case of ^®A1, only a small amount 
of processing occurs on each passage through the radiative zone). Also plotted in Fig. 13 is the 
evolving luminosity of the SEM(1.5) model for comparison. 

For log Tp = 7.76, the CBP flow yields at most 25% of the star’s total energy generation 
rate during the TP-ACB phase. However, this amount is very close to the energy generation that 
would occur in the same radiative layers without CBP; the difference is 1.3% of the star’s total 
energy output. Our proposal to explain ^®A1/^^A1 ratios by CBP requires log Tp extending only 
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to logTp ~ 7.73; this corresponds to 5% of the total energy generation occurring in the CBP flow 
and depletes about 1.2% of the envelope hydrogen. It follows that the CBP necessary to produce 
^®A1/^'^A1~ 10“^ as found in some grains does not have a major effect on the luminosity or the 
energy generation rate. 


3.2. Extra mixing by diffusion parametrization 


We now wish to compare models in which extra mixing is modelled as diffusive transport with 
those in which it is modelled as a mass flow (the present work; Boothroyd, Sackmann & Wasserburg 
1994, 1995; WBS; Boothroyd & Sackmann 1999; Sackmann & Boothroyd 1999; Messenger 2000). 
The mass-flow models have only two parameters, Tp and M. In the case of diffusive models, 
the mixing is modelled either by post-processing calculations which include time evolution of the 
radiative zone, or by inclusion in a full stellar evolution code. In either case, a term is added to the 
stellar evolution model for diffusion between the convective envelope and a region of the radiative 
zone. The parameters used are a diffusion coefficient Umix and a maximum depth of diffusive mixing 
The parameter 5 Mmw is defined in terms of the deepest mass coordinate penetrated, 
the convective-envelope boundary Mbce, and the base of the H shell Mh (see Fig. 1) as 


5Mrr 


Mmix - Mu 
Mbce — Mb 


(24) 


To obtain adequate nuclear processing, the diffusive penetration must extend into the upper part 
of the H shell. To prevent the catastrophe of interference with the region of maximum energy 
generation, a barrier boundary is assumed across which no diffusion is permitted to occur; this is 
specified as (cf. Denissenkov & Weiss 1996; Denissenkov et al. 1998; Weiss, Denissenkov 

& Charbonnel 2000). The two versions of transport - mass flow and diffusive mixing - are thus 
closely related and have similar effects. The relationships between Tp and and between M 

and Umix may be seen as follows. The parameters Tp and S Mmw are directly related through the 
temperature at the position where the diffusive mixing is stopped. This is fixed by the basic stellar 
model. The relationship between Umix and M may be seen as follows. The length scale I and 
time scale t for mixing must satisfy l‘^/Dt ~ 1, so that Umix ~ ^v. For mass flow, the velocity v 
corresponding to M at radial position r in the star is u = M/[47rr^yo(r)], where p{r) is the mass 
density. This yields 


Dr, 


IM 

2 X 10®cm^/s 


M 


I 


10~^Mq/jt I \3 X lO^^cm J \7g/cm' 


-1 


(25) 

(26) 


We have taken for reference the values of r and p at the H shell and I the thickness of the radiative 
zone, all at Stage 5. This is the basic relationship between Umix and M to an order of magnitude. 
For comparison, the value estimated by Denissenkov & Weiss (1996) is Dmix ~ 10®cm^/s in their 
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study of a low-metallicity RGB star. The relations between and Tp for the evolutionary 

stages that we considered for CBP are shown in Fig. 14. In the parametric diffusion treatments of 
Denissenkov and co-authors, 5M^i^ was chosen to be where the H mass fraction was 5-20% below 
the envelope value before the onset of mixing. We indicate the corresponding points in Fig. 14. 
Their prescription results in similar mixing depths to those that we have examined in this paper. 


3.3. The problem of ^®A1 production and Tp 

There remains the issue of whether the calculated effects of CBP for a star may be reasonably 
carried out by choosing a given Tp and M and using a selected radiative-zone structure from a 
single time (stage) during the TP-AGB. A complete calculation of CBP would require inclusion 
of the detailed evolution with time as the stellar structure evolves. Again let us consider a semi- 
analytic treatment to understand the results. The temperature dependence of the processing is 
unimportant for carbon, nitrogen, and oxygen isotopes for all Tp > 7.60. The results for these 
nuclides are thus almost independent of which stellar structure is chosen. We verified this by 
calculations with the numerical model, utilizing radiative-zone structures from Stages 3, 4, 5 and 
6. This is not the case for ^®A1, as can be seen by examining Fig. 3b. For a given Tp, we see 
that the equilibrium ^^Al/^'^Al for different times along the evolution (Eqn. 6) decreases with 
increasing time. This rapid decrease occurs because as the H shell advances outward, the thickness 
(in mass coordinates) of the radiative zone decreases sharply, while the ratio of temperatures at 
the bottom and top of this region, Tp/TpcE (where TpcE is the temperature at the bottom of the 
convective envelope), remains relatively fixed at about 24. The region in which ^®A1 is produced 
accordingly becomes thinner very rapidly as the star evolves, and the amount of ^®A1 produced 
drops. Another contributing factor is that the density in the radiative zone drops steadily during 
the AGB phase, depressing all reaction rates. The final ^^Al/^'^Al value calculated for a given Tp is 
thus drastically altered by the choice of reference state used to provide the radiative-zone conditions. 
That significant ^®A1 production would occur is not an issue. Rather, the question is how does one 
pick an appropriate combination of stellar reference state and Tp for ^®A1 production. Fig. 15 
shows a graph of steady-state ^®A1/^^A1 ratios in the envelope as a function of log(Tp/rH') for the 
different stellar reference states considered. It can be seen that there is a congruent value of Tp for 
each reference state that will produce essentially the same (within a factor of 3) ^®A1/^'^A1 if the 
value of Tp is scaled with Tp as the star evolves. It follows that the Tp required for ^®A1 production 
is within this framework well-constrained and that a GBP model covering all the isotopes of C, N, 
O, Mg, and A1 can be reliably calculated for a given reference state. The particular value of Tp 
is simply congruent to some effective value of Tp/Tp and is not otherwise fixed. This is in accord 
with the approach taken by Boothroyd & Sackmann (1999) and Messenger (2000) for the RGB, 
but with the congruence made evident by ^®A1 production rates. 
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4. Experimental consequences 

We have constructed above a model that should be a good description of AGB stars undergoing 
cool bottom processing, using only two parameters. We now examine the consequences of the rules 
established by this model for observed isotopic abundances. 


4.1. Circumstellar grains 
4 . 1 . 1 . Regime of C/0< 1 

Fig. 16 exhibits the relationships between and in the envelope after first 

dredge-up (curve AF) and after CBP (thin dashed and dotted curves). The available data on oxide 
grains are also shown. A substantial population of the oxide grains can be seen to lie generally 
close to the curve AF for AGB evolution after first dredge-up for stars of varying mass, starting 
from initial solar isotopic composition. (Data are those summarized in Ghoi et al. 1998, and come 
from Hutcheon et al. 1994, Huss et al. 1994, Nittler et al. 1994, Nittler et al. 1997, Nittler & 
Gowsik 1997 and Nittler et al. 1998; curve AF from Boothroyd & Sackmann 1999). Some points 
lie above the curve AF and are plausibly explained by normal AGB evolution without CBP from 
initial states with elevated (Choi et al. 1998). The trajectories for CBP of and 

17o/16o fQj. a, I. 5 M 0 star of solar composition are given by the line from point D to point B, and 
the bounding curve for high M following Fig. 7a. It is seen that there is a population of data 
that lies on this trend and a concentration of data approaching and near to the equilibrium value 
at point B. Following the rules shown in Fig. 7b we see that data in the region labelled H may 
be explained by stars of up to 1.8Mq initial mass. They might equally be explained by an initial 
state of somewhat lower without CBP. Lower-mass stars with solar oxygen ratios directly 

modified by CBP would yield grains which lie in the region bounded to the left by the curve AJB 
and the line AB, corresponding to a range in M and Tp. However, the data in region G can not 
be explained by AGB evolution from solar initial composition, with or without CBP; these must 
involve sources with very low to supersolar initial values and subsolar values - 

presumably representing different stages of Galactic evolution. The issue of how the abundances of 
the rare isotopes and evolve over time in local material and in an average over the Galaxy 
is complex (see Timmes Sz Clayton 1996; Nittler & Gowsik 1996; Alexander & Nittler 1999). It is 
evident that the wide variety of oxygen isotopic abundances found in circumstellar oxides preserved 
in meteorites requires a diverse set of stellar sources with quite different initial isotopic abundances 
and metallicities. Inspection of Figure 10 shows that for C/0 < 1 in the stellar envelopes, there is 
a wide range of Tp that produces the oxygen compositions shown in Figure 16, allowing a plausible 
explanation of a wide range of ^®A1/^'^A1 abundances. The production of ^®A1 is quite adequate to 
obtain values of ^®A1/^'^A1 up to ~ 6 x 10“^, the governing parameter being 7.68 < logTp < 7.73. 
The values of ^®A1/^'^A1< 10“^ cannot be achieved within the framework of standard AGB evolution 
(BGW99) unless there are cases in which essentially negligible ^®A1 production in the H shell is 
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possible and the stellar lifetime is still short enough to contribute to the solar-system grains. 


4-1-2. Regime of C/0> 1 

The regime where carbon or carbide grains may form under equilibrium conditions requires 
C/0 > 1. There are abundant data on and ^^Al/^'^Al in circumstellar carbides recovered 

from meteorites (e.g., Amari, Lewis &: Anders 1994; Hoppe et al. 1994; Hoppe & Ott 1997; Huss, 
Hutcheon & Wasserburg 1997). Most recently the results have been summarized and discussed 
by Amari et al. (2001a,b). In addition to the and ^®A1/^'^A1 data, there are results on 

as well as on Ti and Si isotopes. In the following discussion, we will focus on the so-called 
“mainstream” (MS) grains, “A and B” grains and “Y” grains. These classifications are dominantly 
based on the ratios and (see Amari et al. 2001b for more extensive discussion). 

The mainstream grains are overwhelmingly the most abundant of the circumstellar SiC grains 
found. Figure 17a shows the distribution of grains in the diagram. We first 

consider the mainstream grains. As noted by many previous workers, the of mainstream 

(MS) grains lies in the region expected for the envelopes of AGB stars of approximately solar 
metallicity. The spread of extends from slightly above the solar value of 89 down to 20, 

with most grains between 40 and 80. We also show the whole region (labelled P and Q) that is 
accessible by CBP in the numerical model with C/0 > 1. 

The observed ^®A1/^^A1 versus ^^C/^^C is shown in Fig. 18. The horizontal segment P with 
low ^^C/^^C corresponds to the high-^^N/^^N region also labelled P in Fig. 17a. This is also 
the region in Fig. 11 of low Tp and C/0> 1, running from the solid triangle (no CBP) down to 
12C/13C 

r\j 5 at M > 10 Mq/yi. There is no ^®A1 produced as Tp is low along this segment of 
the curve. The value of ^^Al/^'^Al ~ 8 x 10“^ in region P of Fig. 18 corresponds to the total ^®A1 
produced in the conventional treatment of third dredge-up (BGW99; see Fig. 3a). If we include an 
earlier episode of CBP (prior to TDU) sufficient to give ^®A1/^'^A1 ~ 5 x 10“^, then the envelope 
that is accessible for C/0 > 1 is region R of Fig. 18. The observed ^®A1/^'^A1 for the range of 
12 c/13c found in the MS grains is well represented by CBP. However, the samples with ^®A1/^^A1 
< 6 X 10“^ cannot be explained by the model (independent of C/0). A substantial fraction of these 
grains (and of oxide grains) thus has less ^®A1 than predicted for conventional AGB evolution in 
stars of ~ IMq (BGW99; Mowlavi & Meynet 2000). It is possible that this reflects the normal 
AGB evolution at somewhat lower stellar masses, or that there is some variability in the non-CBP 
^®A1 production in AGB stars as discussed in Sec. 1. 

A fundamental problem with the SiG grains is exhibited by the data shown in Fig. 

17a. Here the majority of MS grains lie below = 10^, and the lower values are near the 

solar value. This is exceedingly difficult to explain with AGB stars in any evolution model that 
assumes solar abundances as the initial state. First dredge-up will destroy roughly half of the 
in the envelope and double the so that the resulting value of should be rsj 10^. Then 

any GBP will only destroy more and produce more increasing the ratio. To explain the 
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high it has been proposed by Huss, Hutcheon & Wasserburg (1997) that the a)^^N 

reaction rates are far greater than the accepted values. If this is the case, then all of the 
calculations presented here are in error. The nuclear data suffer from poorly-known parameters 
for a low-energy resonance, but the most recent work (Champagne &: Pitt 1986) contradicts earlier 
indications of an additional large subthreshold resonance that would produce the much higher rate 
advocated by Huss et al. (1997). If we assume that the calculations presented here are correct 
and that the MS grains represent RGB and AGB evolution, then it is necessary to provide an 
alternative hypothesis. If the MS grains with < 10^ were the products of evolution on the 

RGB prior to first dredge-up then to produce grains with no change from the initial N composition, 
this would require RGB stars that have C/0 > 1 and have done nothing to their N isotopes. No 
evidence for such stars, either theoretical or observational, has been found. If we consider that 
the stars which produced the grains with low were the products of normal first dredge-up 

(DU) and subsequent AGB evolution, then it is necessary that the stars’ intial composition have 
far less than the solar value. In all cases, first DU increases by a factor of ~ 3.7. 

If we assume that the typical initial state had ~ 70, then first DU would give the solar 

value and subsequent GBP would essentially fill the region of the mainstream grains. 

Following these considerations, we carried out calculations in which we increased the assumed 
initial inventory so that = 270 (the solar value) when the star arrives on the AGB. We 

obtained results for all the pertinent isotopes with this initial composition in the numerical model. 
The results for O, C, and Al are of course not affected. The results for and ^^C/^^C are 

shown in Fig. I7b. It can be seen that the whole mainstream region (not just the upper part) is 
covered with this shift in initial We note that adjusting the initial by decreasing 

the amount of while holding all else fixed cannot provide the same adjustment of the final 
This is because first dredge-up with initially solar abundance ratios among the metals 
doubles the mass fraction in the envelope. Adjustments to the initial inventory do not 
affect the gained at first dredge-up, so they can reduce on the AGB by no more than 

a factor of two below the standard models. The SiG grains require a reduction by almost a factor 
of 4. 

It is not clear that the solar value is in any way the proper choice for the ISM sampled 

in forming the stars which are the source of the MS grains. The data on oxide grains clearly 
demonstrate that the stellar sources must come from different molecular clouds which represent a 
wide range in degree of chemical evolution. The requirement that the bulk solar values represent 
a mix of diverse sources is found in many examples. This is particularly the case for rare isotopes 
such as ^^G, and We recall that ^'^0/^*0 in the ISM (Penzias 1981) has been found 

to be distincly different from the solar value (factor of 1.7 ± O.I). We consider that there has been 
evolution of these isotopic abundances over Galactic history, that the abundances of some isotopes 
are quite different today than at the time the solar nebula formed, and that the abundances, 
particularly of very rare or low-abundance isotopes will be distinctly different in different places, 
even at the same time. As is the case for oxygen, the values in the mainstream grains 
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are clear indications that the initial values for some sources cannot be the solar value but must be 
greatly enhanced in This conclusion is unavoidable if the grains are produced by AGB stars. 

Many grains assigned to the MS group exhibit clear s-process signatures in several heavy 
elements which are attributed to nucleosynthesis in AGB stars as noted by Hoppe & Ott (1997) 
and Gallino, Busso, & Lugaro (1997). The more abundant isotopic data on Si and Ti in SiG grains 
show that these elements are dominated by effects associated with galactic evolution (Gallino et 
al. 1994) with small modifications attributed to s-processing. The Si and Ti data are not useful 
for our arguments because they do not strongly reflect the s process. However, the important new 
generation of measurements on individual mainstream SiC grains for many elements show clear 
s-process signatures in several isotopes of different elements (Nicolussi et al. 1997; Nicolussi et al. 
1998; Pellin et al. 2000). In summary, there are compelling reasons to assign MS grains to AGB 
stars, as recognized by many workers. The largest problem with this assignment appears to be the 
data. We conclude that this can most simply be resolved by assuming that the initial 
typically available in molecular clouds ~ 4.6 x 10® yr ago must have ranged from 
~ 70 upward. The “solar” value is considered to be the result of a sampling of material that had 
previously undergone some prior processing in AGB stars. 

The N isotopic compositions determined for the ISM seem to support the idea that 
increases from low values in regions of low astration to higher values in regions of high astration. 
Chin et al. 1999 find ~ 100 for the Large Magellanic Cloud, while higher values of 270 in 

the solar system, 200-600 in the Galactic disk, and > 600 at the Galactic center are found in places 
with presumably increasing astration (see Wilson & Rood 1994 and references therein). It is thus 
plausible that was typically much lower at earlier times and certain that it is variable. 

On the other hand, it is generally claimed (e.g. Wilson & Rood 1994) that there is a gradient of 
14 n/15n with galactocentric distance in the Galactic disk, in a direction suggesting that 
decreases with astration. 

Inspection of Figs. 17a,b and 18 shows the data on Y grains, which are classified as distinctive 
almost solely because their is greater than the solar value. The distributions of these 

grains in both vs. ^^C/^^C and ^®A1/^’^A1 vs. are otherwise the same as those 

for the MS grains. If one considers ^^C dredge-up in models of AGB stars, it is evident that an 
increase in the amount of dredge-up by a factor of 2 or 3 for a small proportion of AGB stars would 
yield a carbon star with Y-grain-like composition. In this case, the Y grains are not distinguishable 
from the more general population of grains produced by AGB stars but reflect some variation in 
C dredge-up. Our preferred interpretation is that they simply reflect the high-end tail of the C/0 
distribution resulting from dredge-up, which may or may not result from variation of metallicity 
(i.e., variation of intial O). The extent to which the variability of dredge-up is possible within 
the framework of AGB models remains to be tested. Other proposals to explain the Y grains are 
based on the dependence of dredge-up on metallicity, as argued by, e.g., Amari et al. (2001b) and 
references therein on the basis of anomalies in the Si isotopic composition. 
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We now consider the A and B grains. As discussed by Amari et al. (2001b), these grains are 
defined by low 10. They are also found to exhibit ranging over more than two 

orders of magnitude and extending to below the solar value. We consider them here as representing 
a group of grains produced by a generic process. It is evident that AGB stars can produce 
as low as 4 at low Tp but they can not produce a wide range in and ^®A1/^’^A1 and 

simultaneously have < 12 and C/0 > 1. It follows that an AGB source at approximately 

solar metallicity can not be the origin of such grains. Further, we do not consider that this condition 
can be directly related to low metallicity. It should also be noted that a considerable number of 
these grains have < 3, ranging down to = 1.57 ± 0.06 and 1.8 ± 0.01 (Amari et 

ah, 2001b). From the existing knowledge of reaction rates, values below 3 can not be achieved on 
either the RGB or the AGB with or without GBP. It follows that either there is some significant 
(x2) error in the reaction rates or that these grains were produced under conditions quite different 
than the CNO burning outlined in Section 2.3 and as discussed by other workers. A considerable 
number of the grains have extending down to = 40, far below the solar value. 

It was suggested by Amari et al. (2001b) that this may be produced by hot H-burning (far above 
AGB shell-burning temperatures), possibly associated with novae. While it is reported that many 
CO novae condense carbon dust, it should be noted that only three of the 31 recent novae listed 
by Gehrz et al. (1998) have been observed to have C/0 > 1 . The formation of different types of 
grains (e.g., carbides or silicates) in the same nova may reflect sequential condensation (see Gehrz 
et ah). The calculated yields for thermonuclear runaways (TNR) in novae almost exclusively have 
C/0< 1 (Starrfield et al. 2001, and S. Starrfield 2001, personal communication). The most recent 
calculations (S. Starrfield 2001, personal communication and Starrfield et ah, in preparation) show 
that C/0 < 1 for white dwarf masses of 1.25 and 1.35 Mq. An extremely high ratio 

(> 1) is predicted in all cases. A search for implanted oxygen with such extreme enrichments 
in graphite or SiC grains would therefore be very important. It is known that extreme values 
of ^^C/^^C and (< 1.5), as well as substantial ^®A1 are predicted for novae yields (cf. 

Starrfield et ah, 2001). Nonetheless, it is difficult to connect the A and B grains as a group to 
novae sources. There are only two reported cases of high ^®A1/^^A1 in SiC grains with extremely 
low ^^C/^^C and They belong to the very small class (6 known) of carbon-rich grains 

attributed to novae on the basis of C, N, and Si isotopic compositions, and they are not A and B 
grains (Amari et ah 2001c). No models of novae give ^^N/^^N far above unity. The yields of C, N 
and O isotopes are all very roughly commensurate (factor-of-ten) so that the wide spread in ^^N/^^N 
cannot be achieved with low ^^C/^^C. Thus, the high ^^N/^^N values found for some A and B grains 
cannot be attributed to novae. The A and B grains appear to suggest some connection with novae; 
however, this cannot be due to the TNR on ONeMg white dwarf models for which nucleosynthesis 
is usually computed. One possibility for a direct nova connection would be if the fragile ^^N could 
be destroyed during the later stages of the explosion. The models of novae so far considered have 
mostly high accretion rates, and it is conceivable that low accretion rates which do not initiate 
full TNR could produce partially processed material with isotopic characteristics resembling A 
and B grains. The true novae represent extremely energetic bursts punctuating a slower, ongoing 
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accretion process. If we assume that the A and B grains were produced by a mixture of accreting 
AGB matter that is not processed in TNR and a component that was produced through TNR, 
then we immediately face the problem that the nova yield ratio Y^/Y^ is typically ~ 0 . 1 , while 
the AGB matter with which it is mixed has almost no A mixing trajectory between AGB 
compositions and nova compositions would therefore move rapidly from the AGB values to those 
with very low but at high ^^G/^^G (see Fig. 17b). Such a mixing curve would not pass 

through the A and B grains. Region P of Fig. 17a,b contains AGB compositions with extensive 
GBP that could mix with a nova composition to produce a mixing curve passing through the A and 
B grains, but only for a small corner of the parameter space where GBP produces ^^G/^^G~ 10 
and 2 x 10^. Since this composition corresponds to low Tp, any ^®A1 found in the grains 

would have to come from the novae. In summary, we have no self-consistent mechanism to propose 
for A and B grains. They appear to reflect different mechanisms than are currently available in 
stellar models. If we assume the classification of A and B grains into a single group to be invalid 
and consider instead subsets of the A and B grains, the problem is not removed. 


4.2. Astronomical implications 

We now consider what evidence concerning cool bottom processing on the AGB should be 
accessible to astronomical observations. As in the case of the pre-solar grains, the potential for 
converting large amounts of G to N is crucial, and we divide up the results into those expected for 
AGB stars with G/0> 1 and with G/0< 1. This discussion applies to stars that have been on the 
AGB long enough to have become G stars without GBP; for the SEM(1.5) model, this corresponds 
to the star exceeding a luminosity of 1.2 x IO^Lq. 

Fig. 10 shows that for almost any value of log Tp > 7.60, M > 10“®MQ/yr can prevent 
carbon-star formation. Low values of M are therefore required if GBP occurs in G stars. Low M 
requires oxygen isotopes with little or no change from the values set by first dredge-up; this means 
i 8 o/i 60 r^j 0.0015 and wide variation of with stellar mass (along the thick-dashed curve 

of Fig. 7). It may further be expected that 10“^ < < 6 x 10“^ and that incomplete 

conversion of G to N in the envelope produces 0.5 < N/0 < 0.8 (cf. Figs. 12a,b) as well as 
12 g/13g> 20. In principle, it is possible to obtain a carbon star with ^^G/^^G< 20 if log Tp < 7.6 
and M is high, but this case constitutes a very small part of the parameter space. Any level of 
^®A1 enrichment is compatible with a G-rich star, depending on Tp. We thus find that contrary to 
usual expectations (e.g. Arnould and Prantzos 1999; Mowlavi & Meynet 2000), normal low-mass C 
stars with low M could produce abundant ^®A1 and they are reasonable targets for gamma-ray line 
observations by high-energy satellite experiments like INTEGRAL. However, low ^^G/^^G requires 
low Tp that is not compatible with significant ^®A1/^'^A1. 

If C/0< 1 in the envelope of an evolved TP-AGB star or planetary nebula, it requires M > 
lO“®M 0 /yr. All of the signs of GBP that become more pronounced with greater M should then 
be present; 8 x 10“^ (Pig. 6 ), 2 x 10"^ (Fig. 12a), and 5 < < 70 
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(Fig. 11). The large quantity of carbon destroyed also requires N/0> 1.5 (Fig. 12b). All of this 
processing should simultaneously destroy ^He. Individual objects may be found with N/0> 3 if 
Tp is very high. Again, since ^®A1 production is independent of M, it should occur equally often 
in 0-rich as in C-rich environments if CBP is active. 

Large enrichments of N/O in planetary nebulae may be important signs of cool bottom pro¬ 
cessing on the AGB. High N/0 is frequently observed and is a defining characteristic of the Type I 
planetary nebulae (cf. Corradi & Schwarz 1995) that make up some 20-25% of observed planetary 
nebulae. Past expectations have been that N enrichment requires HBB (Peimbert &: Serrano 1980; 
van den Hoek & Groenewegen 1997) and thus intermediate-mass stars. CBP would permit N-rich 
planetary nebulae to be made from low-mass (< 3Mq) stars. A search of planetary nebulae for 
other signs of CBP at high M (those associated with 0-rich AGB stars), combined with good stel¬ 
lar mass estimates, would indicate whether the N originated in CBP. Measurements of 
N/0, ^®0/^®0, ^^0/^®0, and ^^C/^^C would provide clear tests of CBP. Evidence for moderate 
N enrichment in a single low-mass planetary nebula has already been found by Pequignot et al. 
(2000). It would be of great interest if more of these objects are found. Since the same H-burning 
cycles that characterize CBP are also active in HBB, it would be necessary to determine initial 
masses of the observed stars and nebulae in order to distinguish CBP and HBB. 


4.3. Galactic inventories 

The extensive processing of stellar envelopes discussed above may have consequences for the 
overall Galactic inventories of at least two nuclides; ^He and ^®A1. Observations of ^He are difficult, 
and have so far been confined to the local ISM and some planetary nebulae (Gloeckler & Geiss 1998; 
Balser et al. 1998; Bania, Rood & Balser 2002). The results indicate that the Galaxy contains no 
more ^He than one would expect from its primordial inventory plus the result of pre-main-sequence 
deuterium burning. It is presently not understood why so little ^He has apparently been made over 
the lifetime of the Galaxy. Low- and intermediate-mass stars should produce large amounts of ^He 
during their main-sequence evolution (e.g., Sackmann & Boothroyd 1999). First dredge-up should 
bring much of this ^He into the stellar envelopes, where subsequent mass loss allows its escape into 
the interstellar medium. Conventional stellar evolution provides no way to destroy this dredged-up 
^He. It has been noted in earlier work (WBS; Charbonnel 1995) that cool bottom processing on 
the red giant branch could provide a path to destroy ^He, using the reactions ^He(^He, 2p)^He 
and ^He(Q!, 7 )^Be. When we include these reactions in our nuclear-reaction network, we find that 
extensive ^He destruction also occurs in the AGB cool bottom processing scenario presented here. 
The results of CBP for ^He are directly analogous to those for there is a threshold temperature 
above which destruction can become rapid. We find a threshold temperature for thorough ^He 
destruction of log Tp ~ 7.65; factor-of-ten depletion may be reached for log Tp = 7.60. Just as for 

the final abundance of ^He in the stellar envelope then depends on the rate and duration of 
mixing, with duration limited by the stellar lifetime. For both ^He and essentially complete 
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(more than a factor of 10) destruction occurs at M > 3 x lO“®M 0 /yr as long as the threshold 
temperature is exceeded. As a result, low-mass stars that undergo cool bottom processing on the 
AGB may not be net producers of ^He over their lifetimes, and this may help explain the apparently 
small amount of ^He production that has occurred. 

An observational test of whether ^He is destroyed in this way may be obtained by spectroscopic 
measurements of in CO lines from AGB stars or planetary nebulae, since both and 

^He destruction should occur simultaneously. The threshold Tp for destruction is lower than 
that for ^He destruction, so destruction is required if ^He is destroyed. The connection between 
^He and destruction is also important for another reason: ^He destruction by GBP must occur 
in essentially all low-mass stars if GBP is to account for the lack of evolution of the ^He inventory. 
This implies that essentially all low-mass stars should destroy most of the in their envelopes, 
and in turn that the preponderance of oxide grains originating in these stars should have 
< 10“^. Examination of Fig. 16 shows that a large fraction of the sources sampled by the oxide 
grains have higher values, suggesting that there might not be sufficient GBP to destroy the extra 
^He produced in low-mass stars. This problem requires further attention. 


The Galactic inventory of ^®A1 is also of interest. The Galaxy contains ~ 3Mq of ^®A1, as 
inferred from observations of the gamma-ray line from its decay (Mahoney et al. 1984; Knodlseder 
et al. 1996). Since ^®A1 has a mean lifetime of 1 Myr, this requires ongoing nucleosynthesis. The 
sites of this nucleosynthesis remain uncertain (e.g., Knodlseder et al. 1999). If the sources of the 
^®A1 are AGB stars (which neither create nor destroy ^'^Al) and if the rate of injection of ^®A1 
into the ISM matches the rate of decay, then the material being ejected into the ISM must have 
^®A 1 /^'^A 1 satisfying 


26Al ^ Mg/t 2 Q 
27 Al “ X 27 M^R’ 


(27) 


where Mq ~ 3Mq is the Galactic 26Al inventory, X’s are mass fractions, M* is the mass of ejecta 
from a single AGB star, and R is the death rate (stars per year) of AGB stars. If X 27 is solar, then 


26Al _ 5.4 X 10-2Mo/yr 


(28) 


in the ejecta (see Dearborn, Lee &: Wasserburg 1988). It is quite difficult to estimate M^R. Con¬ 
sider two cases, with stellar death rates based first on the local density of giant stars and then on 
an estimate of the local white-dwarf birth rate. In either case, we must extrapolate local stellar 
populations into totals for the Galaxy; since the 26Al gamma-ray line emission seems to be concen¬ 
trated in the Galactic disk, we assume that the Galaxy consists only of the thin-disk component 
described in Binney & Merrifield (1998), and we assume that the dying stars (i.e., those shedding 
their envelopes into the ISM and evolving into white dwarfs) make up the same fraction of stars 
everywhere in the Galaxy as they do locally. Using the local density of red giant stars from Binney 
& Merrifield (1998) and assuming that each of these stars sheds its envelope after a lifetime of 
10 ® years, we find R = 1.3 yr If these stars each eject M* ~ ITIq, then the ejecta must have 
26A1/27a1= 4 X 10“2^ from Eqn. 28. Fig. 3a shows that this level could be produced by GBP 
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with log Tp ~ 7.75. Extensive CBP at high Tp might permit AGB stars to be a significant sonrce 
of the Galactic ^®A1. However, the preponderance of circnmstellar grains appear to be produced 
in environments where ^®A1/^^A1 is much below 4 x 10“^. If we use the white-dwarf birth rate 
of Drilling & Schonberner (1985) and we extrapolate the local rate of (2 to 3) x 10“^^ pc“^ 
to a total Galactic rate, we obtain R = 0.04 to 0.06 yr“^. This is far too low to contribute to 
the Galactic inventory. We conclude that it is unlikely that a significant fraction of the ^®A1 ob¬ 
served in gamma-ray line emission in the Galaxy was made in AGB stars. This conclusion has 
also been reached on the basis of the apparently clumpy distrubution of A1 in the Galaxy (Diehl 
and Timmes 1998). However, cool bottom processing in AGB stars may provide sufficient ^®A1 to 
permit detection in selected AGB stars with very low and high N/0. 


5. Conclusion 

The work presented in this paper demonstrates that extensive predictions of the possible range 
of envelope compositions and isotopic abundances of AGB stars may be made in a model of cool 
bottom processing parametrized only by a mass-flow rate M and a maximum temperature Tp using 
a simple calculation. 

This calculation uses results from a stellar evolution model calculated without GBP to provide 
the basic conditions needed to determine the results of additional processing by extra mixing. It 
appears that significant composition changes have negligible feedback on the stellar luminosity. 
As a result, the post-processing approach adopted here - in which CBP is calculated separately 
from stellar evolution - appears to be reasonable. Post-processing calculations allow the range 
of possible effects of extra mixing to be determined easily by exploring the space of Tp and M 
values once stellar structures have been determined for an arbitrary star. The results reported in 
this work are based on a grid of 231 pairs of the two parameters. This covers the entire range 
that is both plausible and potentially interesting, using an unsophisticated code and about one 
day of processing time on a personal computer with 930 MHz clock speed. Moreover, given only 
the stellar structure and reaction rates, one can obtain almost quantitatively correct analytical 
approximations to these results as functions of M and Tp for C/0 (Eqn. 22 and Eig. 10), ^®A1/^'^A1 
(Eqn. 6 and Fig. 3a), (Eqn. 13 and Fig. 4), and (Eqn. 16 and Fig. 5). Given 

the uncertainties concerning the exact mass, evolutionary state, mass-loss rate, and extra-mixing 
parameters corresponding to the star that produced an observed presolar grain or planetary nebula, 
the analytic approximations are reasonably accurate and display in a clear way the rules that govern 
the outcome of CBP. 

We found also that given the uncertainties (particularly the unknown mixing mechanism, which 
will set the mixing depth), the detailed evolution of the stellar structure below the stellar envelope 
(e.g., the radiative zone) is relatively unimportant. The results from a complete stellar evolution 
code without extra mixing are necessary to specify the timescales, the evolution of envelope mass, 
the action of dredge-up, and a representative structure for the region in which nuclear processing 
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takes place in the CBP model. However, the detailed evolution of this region need not be followed 
to obtain useful results. The major effect of the evolution of this region is to change the Tp that 
corresponds to a given amount of processing; holding both the structure of this region and Tp fixed 
is equivalent to holding Tp/Tn fixed in a model that follows the detailed evolution, provided that 
both follow evolution of the envelope in the same way. 

The effects of CBP on stars of different masses or metallicities (particularly low metallicity) 
can readily be determined once a conventional non-CBP stellar model and radiative-zone structure 
are obtained from conventional codes. 
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Chieffi and M. Limongi, as well as extensive comments and suggestions from J. Lattanzio and an 
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to thank the Lunatic Asylum and Caltech for hospitality. This work was supported by NASA 
grant NAG5-11725 and by the Italian MURST contract Cofin2000. Caltech division contribution 
8769(1080). 


A. Reaction Rate Uncertainties 

Most of the reaction rates used in this paper are known to within 5-40%, according to the 
NACRE evaluation (Angulo et al. 1999). There are a few cases of larger uncertainty, which we now 
examine. In cases where more or less reliable error estimates exist, the NACRE upper and lower 
limits to the reaction rates are often found by multiplying by two the size of the statistical error 
bars one derives directly from the experimental reports. We consider the rates at a temperature of 
5 X lO'^ K as representative of the crucial temperature range. 

Eirst, the NACRE compilation quotes uncertainties of a factor of two for the rate of ^®0(p, a)^^N. 
This affects only the rate of destruction at low Tp (< 10^'®^ K). Above this temperature, in¬ 
creasing the destruction rate changes nothing because even with the recommended rates, a single 
pass through the radiative zone suffices to greatly reduce the abundance in the circulating 
matter. Likewise, decreasing this rate will have no effect as long as the lower rate still suffices to 
destroy all in a single pass; the threshold Tp for this condition is also about lO'^'®^ K. Below 
log Tp = 7.62, factor-of-two changes in the 0-destruction rate translate into factor-of-two changes 
in the final 

The most important uncertainties are those concerning reactions which produce and destroy 
^®A1. The NAGRE compilation quotes upper and lower limits for the rate of ^®AR(p, 7 )^'^Si which 
are a factor of 47 higher and a factor of 86 lower, respectively, than the recommended rate. It should 
be kept in mind that the NACRE recommended rate at these temperatures was determined mainly 
by setting resonance strengths to be 1/10 of upper limits derived from the shell-model calculation 
of Champagne, Brown & Sherr (1993; these are considerably lower than the experimental upper 
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limits). Changing from the NACRE recommended rate to the NACRE upper-limit rate reduces 
the equilibrium ^®Al/^^Mg in circulating material from near 100 to around unity. As a result, low 
M then allows this ratio to equilibrate in the circulating material, so that the proportionality of 
^®A1 production to M~^ in Eqn. 4 breaks down, and the amount of ^®A1 returning to the surface 
becomes both lower and independent of M. Although the final ^®A1/^^A1 in the envelope is smaller 
in this case, it still remains high enough to account for all oxide grains with only a very slight 
increase in the maximum Tp. 

The rate of the ^®A1 production reaction, ^^Mg(p, 7 )^®AR, is quoted with an upper limit 75% 
above the recommended rate and a lower limit 41% below the recommended rate. The NACRE 
rate and uncertainties for this reaction are based on resonance strengths extracted by Iliadis et ah 
(1996) from pre-existing data using improved DWBA methods. Except at M < lO“'^M 0 /yr, where 
the independence of the envelope ^®A1 enrichment rate from M breaks down, the result of varying 
the ^®Mg(p, 7 )^®A 1 rate within the NACRE limits is to scale the final ^®A1/^^A1 downward by as 
much as 40% or upward by up to 75%. 

The only other large uncertainties (greater than a factor of two) in the reaction rates included 
in our code are for the ^'^Al-producing and -destroying reactions ^®Mg(p, 7 )^^A 1 , ^^Al(p, Q;)^^Mg, 
and ^'^Al(p, 7 )^®Si. The first of these reactions has an upper limit in the NACRE compilation which 
is 15 times the recommended rate (derived mainly from experimental work by Champagne et ah 
1990) and a lower limit which is 1/27 the recommended rate. The upper and lower limits to the 
{p, a) destruction rate differ from the recommended rate by factors of 10 and 5, respectively. To 
examine the allowed effect on ^®A 1 /^^A 1 due to ^^Al production and destruction (which is negligible 
when using the recommended rates), we ran models with all of the ^^Al-producing reaction rates 
(including the ^®Al-destroying ^®Al(p, 7 )^^Si(/ 3 '''i/)^^Al) set to the NACRE upper limits and the 
^^Al-destroying reaction rates set to their NACRE lower limits. This made virtually no difference 
to the final ^®A1/^^A1 beyond the effects of ^®Al(p, 7 )^'^Si on ^®A1 production described above. The 
corresponding calculation with the limits switched to minimize the final ^^Al abundance showed 
no significant difference from that using the recommended rates. 
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Fig. 1.— Schematic diagram of the cool bottom processing model. Material taken from the 
envelope circulates slowly down into the radiative zone, where it undergoes nuclear processing 
at temperatures near Tp. It then returns to the envelope, where it is rapidly mixed with the other 
envelope material. Relative thicknesses of the various layers do not reflect physical dimensions. 
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Fig. 2.— Evolution of the boundary of the C/0 core Mcore, location of the H shell Mh, and 
convective-envelope boundary Mbce for the 1.5Mq solar-metallicity model of Straniero et ah 
(1997) [SEM(1.5)], as functions of time since the He flash. Vertical lines indicate times selected to 
provide radiative-zone structures for our CBP calculations (cf. Table 1). 
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Fig. 3.— a) Final ^^Al/^^Al as functions of Tp from the numerical model running for the 8.5 x 10^ 
yr period during TDU (Stage 5 to Stage 6), compared with the corresponding semi-analytic result 
for the CBP steady-state abundance (curve 5 of panel b). “Extreme oxide grains” indicates the 
highest value found for a meteoritic oxide grain, b) Envelope ^^Al/^^Al at steady state ( Eqn. 7). 
The production rate P changes as the star evolves; each curve shows the value of for a 

particular choice of radiative zone structure taken from the stellar evolution model. Curves are 
labelled according to stage numbers in Table 1. The end of each curve indicates Tp at that time. 
The “std evolution” line indicates the level found for TDU production of ^®A1 by BGW99. 
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Fig. 4.— Fraction of surviving in the final stellar envelope, as a function of M, after 8.5 x 10^ 
yr. The “analytic” curve is based on a stellar envelope of O.SMq and complete destruction of 
in the circulating material (Eqn. 13) . The solid curves correspond to processing for the same 
amount of time in the numerical model, at four different values of Tp. 
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Fig. 5.— ratio in the final stellar envelope after processing throughout the TDU period. 

The analytic curve (Eqn. 16) reflects an initial envelope with composition set by first dredge-up, 
mixing with processed material with the composition 0.0011. The other curves are 

results from the corresponding numerical model, at 4 different values of Tp. 
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Fig. 6.— Final envelope compositions from the numerical model with the Stage 5 radiative-zone 
structure. The grid consists of curves of constant M (dotted lines) and constant Tp (dashed lines). 
Superimposed over the grid are the oxide grain data summarized in Choi et al. (1998); filled circles 
are those with ^®A1/^^A1> 2 x 10“^. In principle, one can read values of M and Tp off of this and 
subsequent figures; for example, a grain with ^®A1/^^A1= 3.5 x 10“^ and 0.00087 (at the 

square) corresponds to the final product of CBP with log Tp = 7.728 and \og{M/M q) = —6.4. 
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Fig. 7.— Final envelope compositions of a I.SMq star after CBP (10“®MQ/yr < M < lO“^M 0 /yr, 
7.44 < logTp < 7.76, notations as in Fig. 6). a) Point D is the composition after first dredge-up. 
The triangle B is the composition reached if CBP continues for a long time at log Tp > 7.62. 
The calculated compositions all lie near a two-component mixing curve connecting points D and 
B. Compositions to the left of this curve indicate incomplete processing. The thick dashed curve 
ACDEF indicates the results of dredge-up during the red giant phase of stellar evolution for various 
initial stellar masses (cf. WBS). b) trajectories resulting from CBP for two arbitrary stars with 
O compositions of A' and A" after first dredge-up. Mixing lines point to the equilibrium 
value. Range of M and Tp is the same as in panel (a). 
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Fig. 8.— Equilibrium abundance ratio as a function of temperature, which depends only 

on the reaction rates. The high values on the left side of the graph are inaccessible in CBP because 
the reactions are too slow to establish equilibrium at logTp < 7.6 (See Table 1). 
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t (Myr since start of TDU) 

Fig. 9.— a) Time evolution of C/0 in SEM(1.5) (no CBP). Also shown is the envelope mass as 
a function of time. The C/0 curve labelled “with CBP” is the evolution of C/0 with time for a 
star with CBP. The N/0 curve is for the same CBP model. Note that CBP adds to N/0 in the 
envelope at the expense of C/0, b) Evolution of ^^C/^^C and ^®A1/^^A1 with time in the numerical 
model of concurrent CBP and TDU, at two choices of M and compared with the no-CBP model. 
Arrows indicate the general direction of evolution. The jumps in composition are due to dredge-up 
episodes, and motion toward lower ^®A1 reflects radioactive decay. Increasing Tp moves the end 
points of these curves to the right. Increasing M moves the end points downward. 
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Fig. 11.— Final compositions of the envelope where CBP stops immediately before the last thermal 
pulse, after which the last dredge-up episode adds Tp increases clockwise along the constant-M 
curves; dots on the M = lO“®M 0 /yr curve indicates log Tp values in steps of 0.016, and points 
with logTpi = 7.616, logTpi = 7.632, and logTps = 7.648 are circled and labelled. The solid 
triangle indicates the final composition of a 1.5Mq star without CBP. Shaded region contains 
M < lO“®M 0 /yr, approximately the region of carbon star formation, while the large open triangle 
encloses compositions which have C/0> 1 and ^®A1/^'^A1> 10“^. Note that low-^^C/^^C carbon 
stars are restricted to a small region in the lower right-hand corner, near the solid square. EQ is 
the limiting equilibrium value. Bounding curve “EQ -|- final TP” corresponds to added by the 
final dredge-up episode, mixing with an envelope of composition EQ. 
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Fig. 12.— Final envelope compositions for the numerical model, a) decreases from an 

initial value of 6.3 x 10“® toward a final value ~ 3 x 10“® for all Tp, but with an effectiveness 
that varies mainly with M. The arrow at the top indicates the envelope composition after first 
dredge-up. b) final values of N/0 and C/0. At most temperatures, this reflects varying degrees of 
conversion of C into N; the two curves with the log Tp > 7 reflect additional destruction of O to 
make more C and N. The point labelled “solar” indicates solar-system composition. 
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Fig. 13.— Solid curve; mean luminosity in the CBP circulating regions as a function of Tp, 
estimated as described in the text (read bottom axis). Dashed curve; luminosity of the SEM(1.5) 
model without CBP, as a function of time (read top axis). 
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Fig. 14.— Relation between the parameters (5Mmix and Tp, at each of the 6 selected evolutionary 
stages of the radiative zone (cf. Table 1 and Fig. 3b). Along the Stage 5 curve, the thick broken line 
indicates Tp considered in this study; the connected black dots indicate depths chosen to have H 
depleted by 5-20% as assumed for S Mmw in Weiss, Denissenkov & Charbonnel (2000). This shows 
the essential equivalence of the mass flow approach [M,Tp) and the diffusive mixing approach. 
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Fig. 15.— Steady-state {t 2 qP/Y^) calculated for the envelope assuming the whole 

evolution takes place with the radiative-zone structure fixed at Stages 1-6, as labelled. The results 
are as in Fig. 3b but in terms of log (Tp/Tn)- The ^®A1 production is governed by the last stages. 
Note that for Stages 4, 5, and 6, abundant ^®A1 is produced for \og{Tp/Tp) > —0.10 and that 
PocM-^ 
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Fig. 16.— Final envelope and with oxide-grain data (summarized in Choi et ah 

1998) superimposed. The dotted curve AJB indicates the leftmost limit of compositions available 
starting from solar initial composition without dredge-up. Region G is therefore not accessible by 
CBP for a star with solar initial composition. 
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Fig. 17.— and of SiC grains compiled by Amari et al. (2001b). Symbols are: o 

A+B grains, □ mainstream (MS) grains, o Y grains. The histogram at the bottom is the frequency 
of each type (arbitrary linear units), a) Regions P and Q contain all compositions with C/0> 1 
produced during TDU and CBP; Q is with ^®A1/^'^A1> 10“^. b) same model as (a), but for 
14 n/15n = 270 when the star arrived on the AGB, corresponding to an intial ~ 70 when 

the star formed. The dot-dashed curve shows compositions that result when a typical low-^^N/^^N 
composition from an AGB star mixes with ejecta from a nova outburst. Points along the curve 
indicate the fraction of the composition made up of nova ejecta (nova composition from Starrfield 
et al. 2001) 
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Fig. 18.— C and A1 isotopic compositions of SiC grains compiled in Amari et al. (2001b; symbols 
as in Fig. 17). The region P-Q contains all compositions produced with TDU and CBP. For all 
values of Tp and M considered (except for the small part of region P with 50, these 

are the same regions as in Fig. 17). Region R is analogous to P-Q, but for models which assume 
that extensive CBP produced ^®A1/^^A1= 10“^ and ^^C/^^C=3 before TDU, and then continued 
throughout TDU. 




